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ABSTRACT 

We use high resolution direct numerical simulations (DNS) to show that helical tur¬ 
bulence can generate significant large-scale fields even in the presence of strong small- 
scale dynamo action. During the kinematic stage, the unified large/small-scale dynamo 
grows fields with a sh ape-invariant eigen fu nction , with mo st power peaked at small 
scales or large k, as in ISubramanian fe Brandenburg) (I2Q14T) . Nevertheless, the large- 
scale field can be clearly detected as an excess power at small k in the negatively 
polarized component of the energy spectrum for a forcing with positively polarized 
waves. Its strength B, relative to the total rms field B rms , decreases with increasing 
magnetic Reynolds number, Rejvi- However, as the Lorentz force becomes important, 
the field generated by the unified dynamo orders itself by saturating on successively 
larger scales. The magnetic integral scale for the positively polarized waves, charac¬ 
terizing the small-scale field, increases significantly from the kinematic stage to sat¬ 
uration. This implies that the small-scale field becomes as coherent _as possible for 
a given forcing scale, which averts the ReM-dependent quenching of B/B rms . These 
results are obtainedfor 1024 3 DNS with magnetic Prandtl numbers of PrM = 0.1 and 
10. For PrM = 0.1, B/B lms grows from about 0.04 to about 0.4 at saturation, aided in 
the final stages by helicity dissipation. For PrM = 10, B/B rms grows from much less 
than 0.01 to values of the order the 0.2. Our results confirm that there is a unified 
large/small-scale dynamo in helical turbulence. 
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1 INTRODUCTION 

Astrophysical systems like stars and galaxies host magnetic 
fields that are coherent on the scale of the system itself. 
These fields are thought to arise due to the action of a tur¬ 
bulent dynamo, whereby helical turbulence combined with 
shear amplifies and maintains fields coherent on scales larger 
than the scales of random stirring. Indeed, the scales of the 
stirring like convective scale in the Sun or the supernova- 
induced turbulent scales in galaxies, are much smaller than 
the coherence scale of the large-scale field. A dynamo which 
amplifies fields on scales larger than that of the stirring is 
referred to as a large-scale or mean-field dynamo. 

There are two major potential difficulties associated 
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with mean-field dynamos. One is that small-scale helical 
fields which are produced during mean-field dynamo ac¬ 
tion, due to magnetic helicity conservation, go to quench 
the dynamo. Thus, they have to be eliminated from the 
dynamo active region by so me form of magnetic helicity 
flux, t o avoid such quenching (IBr andenburg fe S ubramanianl 


l2005l : IBrandenburg et all 1 20121 : Blackman 2015l l. Equally 


important is the fact that, while mean-fields are being 
generated, the same turbulence, for a large enough mag¬ 
netic Reynolds number, ReM, also generically lead to the 
small-scale or fluctuation dynamo. The fluctuation dy¬ 
namo rapidly generates magnetic fields coherent on scales 
of the order of or smaller than the outer scales of the 
turbulence, and in principle, at a rate faster than the 
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to whether and how the mean-field dynamo operates in the 
presence of such rapidly growing magnetic fluctuations. 

This issue was p artially addressed by 
ISubramanian fc Brandenburg! (hoidt ) (hereafter SB14), 
throug h direct numerical simulations (DNS) and by analyz¬ 
ing the iKazantsevI (119681 4 model, focusing on the kinematic 
regime. They showed that in this regime, the magnetic 
energy spectrum grows as an eigenfunction, i.e. at each 
wavenumber fc the spectrum grows with the same growth 
rate. Nevertheless, there is indeed evidence for this large- 
scale field generation in the horizontally averaged large-scale 
field, which can also be seen as excess power at small k 
in one of the oppositely helically polarized components. 
However, they also found that the relative strength of the 
large-scale or mean-field component compared to the rms 
field, in the kinematic stage, decreases with increasing ReM 
like Re^ 3 " 74 for larger values of ReM- From both an analysis 
of the Kazantsev model including helicity jKazantseyl 
196^_ |\4nnsj^in_fc Kichatinovl Il986 : ISubramanianl 1 19991 : 


Boldyrev et al. 20051 ) and the DNS. SB14 showed that this 


is a result of the magnetic energy spectrum peaking on 
small resistive scales even in the presence of helicity. If 
such a feature persisted on saturation, it would be difficult 
to explain the prevalence of large-scale fields. Of course, 
as the dynamo-generated field grows, the Lorentz force 
will become important, first at small scales, and saturate 
the growth of small-scale fields. It is then important to 
determine whether the mean-field dynamo can continue to 
grow large-scale fields as the fluctuation dynamo saturates. 
And can this large-scale field become more dominant at 
saturation, independent of ReM? Our aim in this paper is 
to answer these questions. 

For this purpose we have run DNS of magnetic field 
growth in helically driven turbulence in a periodic domain, 
with resolutions up to 1024 3 mesh points. These simulations 
are designed to adequately capture the dynamics of scales 
both smaller and larger than the forcing scale, and run from 
the kinematic regime to nonlinear saturation. The next sec¬ 
tion presents the simulations that we have carried out to 
use for our analysis. Section 3 sets out the results of our 
analysis to determine the evolution of both the large- and 
small-scale fields generated by helical turbulence. Section 4 
presents a discussion of these results and the last section, 
our conclusions. 


2 SIMULATING LARGE-SCALE DYNAMOS 

To study the growth of the large-scale or mean-field dy¬ 
namo in the presence of a small-scale or fluctuation dynamo, 
we have run a suite of simulations of helically driven tur¬ 
bulence using the Penc il CodeQ jBrandenburg fe Doblerl 
l2002l ; lBrandenburg|2003|) . The continuity, Navier-Stokes and 
induction equations are solved in a Cartesian domain of a 
size (27r) 3 on a cubic grid with N 3 mesh points, adopting 
triply periodic boundary conditions. The fluid is assumed 
to be isothermal, viscous, electrically conducting and mildly 


1 https://github.com/pencil-code 


Table 1. Summary of runs discussed in this paper. Here the val¬ 
ues of rirms are from kinematic phase, whereas B rm s and B refer to 
values from the saturation phase after ~ 624 eddy turnover times 
for Runs A, B, and C and after ~ 463 eddy turnover times for 
Runs D, E, and F; the eddy turn over time is given by l/(u r msfcf). 


Run 

PrM 

Re M 

^rms 

-Rrms 

B 

B/B lms 

N 3 

A 

0.1 

330 

0.135 

0.085 

0.033 

0.38 

1024 3 

B 

0.1 

160 

0.130 

0.092 

0.046 

0.49 

256 3 

C 

0.1 

65 

0.130 

0.092 

0.055 

0.59 

256 3 

D 

10 

3375 

0.135 

0.078 

0.017 

0.22 

1024 3 

E 

10 

1575 

0.126 

0.072 

0.019 

0.27 

256 3 

F 

10 

665 

0.133 

0.082 

0.026 

0.32 

256 3 


compressible. The governing equations are given by, 

j|lnp»-V-u, (1) 

-2-U = Cg Vlnp + - J x B + EvIbc + /, (2) 

Ut p 

d 

— A = ux B - pp 0 J. (3) 


Here p is the density related to the pressure by P = pc 3 , 
where c s is speed of sound. The operator D/Dt = d/dt+u-X7 
is the Lagrangian derivative, where u is fluid velocity field. 
The induction equation is being expressed in terms of the 
vector potential A, so that JB = V x A is the magnetic field, 
J = V x .B/po is the current density and po is the vacuum 
permeability (po = 1 in the DNS). The viscous force is given 
as F v i S c = p _1 V • 2upS, where v is the kinematic viscosity, 
and S is the traceless rate of strain tensor with components 
Sij = ^(ui,j +Uj,i) — V u. Here commas denote partial 
derivatives. The forcing term / = /(x, t) is responsible for 
generating turbulent helical flow. The forcing is maximally 
helical as described in SB14. In Fourier space, this driv¬ 
ing force is transverse to the wavevector k and localized in 
wavenumber space about a wavenumber fcf. It drives vorti¬ 
cal motions in a wavelength range around 2n/k{, which will 
also be the energy carrying scales of the turbulent flow. The 
direction of the wavevector and and its phase are changed 


at every time step in the simulation making the force 5- 

correlated in time; see Brandenburg ( 

2001 

): Haugen et al. 

(20041): Subramanian & Brandenburg 

2014 

) for details. 


For all our simulations, we choose to drive the mo¬ 
tions at a wavenumber fcf = 4. This choice is motivated by 
the fact that we wish to resolve both the small-scale mag¬ 
netic field structures in any turbulent cell and at the same 
time include scales larger than the flow (with k < fcf). The 
strength of the forcing is adjusted so that the rms Mach 
number of the turbulence, u rms in the code (where veloc¬ 
ity is measured in units of the isothermal sound speed), is 
typically about 0.1. This small value implies also that the 
motions are nearly incompressible. We define the magnetic 
and fluid Reynolds number respectively by ReM = Wrms/pfcf 
and Re = Mrms/i'fcf, where p and v are the resistivity and 
viscosity of the fluid. The magnetic Prandtl number is de¬ 
fined as PrM = ReM/Re = v/p. 

For the simulations reported here we take PrM = 0.1 
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and PrM = 10; the former value was used in most of the 
DNS in SB14. This was motivated by the fact that for such 
small values of PrM, the non-helical sm all-scale dynamo wa s 
expected to be much harder to excite (llskakov et al.ll200it) . 
which would then provide a better chance of seeing evidence 
for the large-scale field in the kinematic stage. SB14 however 
found an efficient small-scale dynamo even for PrM = 0.1. 
This is related to the fact that in SB14 the forcing wavenum- 
ber was chosen to b e k{ = 4, while in the earlier work of 
llskakov et al.l (120071 ) it was between 1 and 2. Furthermore, 
the disappearance of the small-scale dynamo for k{ = 1.5 
is related to the bump in the spectrum near the dissipa¬ 
tion waven umber, which is known as the “bottleneck phe¬ 
nomenon” rtFalkovichl 1 1994 1. In the nonlinear regime, how¬ 
ever, this bump is suppressed by the magnetic field and 
therefore the strengt h of the small-scale dynamo is nearly 
independent of PrM (lBrandenburdl201 ll ). As we would like 
to compare with the kinematic results of SB 14 and extend 
it to the nonlinear regime, we consider first the value of 
PrM = 0.1. For both values of PrM, we have studied a range 
of ReM- The fiducial simulation of a helical dynamo con¬ 
sidered in this paper has a resolution of 1024 3 mesh points 
and Rcm = 330 (hence Re = 3300), which we from now on 
refer to as Run A. The other run with a similar resolution 
of 1024 3 mesh points has PrM = 10 and is referred to as 
Run D. A summary of different runs used in this paper is 
given in Table [T| 

A quantity that was also evaluated in the DNS of SB14 
and which helps to understand how the growth rate of the 
fluctuation dynamo changes as compared to the mean-field 
dynamo is the growth rate as a function of k, A (k), of the 
magnetic energy residing in each k defined in the following 
manner. We first extract the evolution of magnetic energy at 
a given k, Mk(t) from Bm)^, t) and perform a running aver¬ 
age with a suitable window, to smooth the evolution curve, 
Mk(t). This window is chosen such that the smoothing re¬ 
duces the noise sufficiently but at the same time does not 
produce any new false features in the curve. Subsequently, 
to determine A, we take the logarithmic derivative given by 


1 d Mfe(t) 

M k {t) At 


(4) 


We determine the A for each k, to obtain A (k). 



Figure 1. Evolution of Eyi(k,t) for Run A (solid black lines) 
together with (dotted blue), with the final curve in solid 

blue. The spectra are at regular intervals of 100 code time units, 
starting at t = 100. 



Figure 2. Evolution of Mfc(t) for k = 1,4. 50, 200, and 500 has 
been shown for Run A. 


3 RESULTS 

We first consider the results from one of our higher resolu¬ 
tion simulations of a turbulent helical dynamo with 1024 3 
mesh points (fiducial Run A) from the kinematic stage up to 
nonlinear saturation. Most of our results have been obtained 
using the spectral data from the DNS. 

3.1 The evolution of spectra 

In Figure [T| we show the evolving kinetic and magnetic 
energy spectra, E\^{k) and Eyilk), of the helically driven 
dynamo, in dotted blue and solid black lines respectively. 
For a maximally helical velocity field, the mean-field a 2 dy¬ 
namo in a periodic box, is expected to grow fields initia lly 
at a w avenumber k — fcf/2 = 2 (see Brandenburg et abl 
(|2002|) and SB 14), which will then move to even larger scales 


(smaller k) in the saturated state dBrandenburg||200ll ). On 
the other hand, the fluctuation dynamo is expected to be 
active at scales equal to and below the forcing scale (larger 
wavenumbers). In the early kinematic stage, a single com¬ 
mon eigenfunction is seen growing in a self similar manner, 
with the magnetic energy spectrum Eyi{k,t) increasing to¬ 
wards the smaller resistive scales. This is similar to what 
is seen in the DNS of large-scale dynamo action by SB14, 
during the kinematic stage. Also the slope of Eyi[k) is close 
to k 7 / 6 , which agrees with a result derived by SB14 using 
Kazantsev model including helicity. 

At late times, when the dynamo saturates, the slope 
of EM(k,t) as a function of k flattens first, with the peak 
of the spectrum shifting secularly to smaller and smaller 
wavenumbers (larger and larger scales). The subsequent sat¬ 
urated spectra develop two peaks, one at the forcing scale 
and another at the largest scale or smallest wavenumber, 
k = fci = 1. 
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Figure 3. Evolution of M/.(i) for k = 1, 4, 10, 50, and 200 has 
been shown for Run D. 



Figure 5. Growth rate A (k) for Run A as a function of k at times 
t = 200, 300, 400, and 500. 




Figure 4. In the top panel, spectra of E M for Run A are in solid 
curves and E ^ are in dotted curves at three times, t = 200, 500, 
and 900. In the bottom panel the corresponding growth rate, A(fc) 
is shown. 


than for Run A (Pi~m = 0.1). This is because for PrM > 1 
(and also large ReM), firstly, the smallest eddies whose 
scale-dependent ReM is greater than the critical value for 
small-scale dynamo action, have a shorter turnover time. 
Secondly, the larger value of ReM for Run D (by an or¬ 
der of magnitude, compared to the PrM =0.1 run), could 
make the dynamo in that run more efficient - even for a 
similar eddy turnover rate. Thus, the helical dynamo in 
Run D turns to saturate at a much earlier time compared 
to Ru n A. As is w ell know n starting from the wor k of 
iBrandenburd (l200lll fsee al solBrandenburg fe S ubramanianl 
( 2005ll : iBrandenburg et all ll2012li : Blackman ( 2015l l for re- 
views), this growth of large-scale field requires small-scale 
helicity to be lost from the system. In the present context 
of a completely homogeneous dynamo with uniform energy 
density of the large-scale field, such a loss is purely due to 
resistivity, whereas more realistically it would be aided by 
helicity fluxes out of the dynamo active region. Nevertheless, 
the magnetic field evolution reflected in Figures [l] and [2] goes 
to show that even after the fluctuation dynamo saturates, 
the mean-field dynamo continues to grow large-scale fields, 
provided also that small-scale magnetic helicity can be lost. 


To make this saturation behavior clearer, we show for 
Run A in Figure [2] the evolution of magnetic energy resid¬ 
ing in different scales. The evolution of Mk is shown for 
wavenumbers k = 1, 4, 50, 200, and 500. It can be seen that 
Mk first grows almost exponentially with similar slopes for 
all k in the kinematic stage, before saturating. For higher 
wavenumbers, Mk stops growing and turns to saturate at 
lower strengths, and earlier than magnetic energy at lower 
wavenumbers. The M 4 and M 50 modes (dotted blue and 
dashed blue lines), turn to saturate at around t = 525 and 
450, respectively, as compared to Mi in solid black, which 
has not saturated even at late times. In fact, the Mi mode, 
which reflects the operation of the large-scale dynamo, is 
seen to be still growing, and has a distinctive positive slope 
compared to the saturated flatter curves at other wavenum¬ 
bers. 

A similar picture is seen also for Run D (PrM = 10), 
as shown in Figure [3] Here, the growth rate is much larger 


An interesting behavior to note is that the Mi mode is 
growing together with all other modes at about the same 
rate and turns off to saturate almost along with M 4 . From 
this a picture emerges where there is one unif ied large- and 
sma ll-scale dynamo in such helical turbulence IlSubramanianl 
11999h . which simply grows fields on all scales together, and 
saturates at successively larger and larger scales (smaller 
and smaller wavenumbers). We return to this idea below. 

3.2 Polarization spectra and wavenumber 
dependent growth rate 

For the large-scale dynamo action which arises in heli¬ 
cal turbulence, the turbulent emf, in a two-scale model, 
is expected to generate oppositely signed small- and large- 
scale magnetic helicities. Thus, one way of distinguishing 
large- and small-scale fields would be to compute spec¬ 
tra from the field split into positively and negatively po- 
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larized components, defined as dBrandenburg et al.l 120021 ; 

iBrandenburg fc Subramanianll2005l~l 

E±(k,t) = 5 [EM(k,t) ± \kHwi{k, t)] . (5) 

This would enable one to see a clearer signature of the large- 
scale field and its evolution from the kinematic stage to non¬ 
linear saturation. In the top panel of Figure [3] we show 
for Run A the spectrum of the two oppositely polarized 
field components (depending on lielicity), at three 

times, t = 200, 500, and 900. The spectra of the negatively 
polarized field, are shown as solid lines and those 

of the positively polarized field are shown as dotted lines. 
Also the continuous ordering of the field, and the continued 
growth of large-scale field even when the small-scale field 
saturates, can also be seen by examining the wavenumber- 
dependent growth rate, A (k). In the bottom panel of Fig- 
ure[4j the growth rates A (k) corresponding to E^(k,t) have 
been plotted. The black curves in both panels are at t = 200, 
in the kinematic stage. We see that at this time A (k) is nearly 
uniform across k in both polarized components, indicating 
in the presence of an eigenfunction which is growing at the 
same rate at all scales, both large and small. (This can also 
be seen from Figure [5]) We see that both polarized spectra 
extend over all k, and in fact overlap at large k, indicating 
that the magnetic field has little helicity on such large k. 
However even at this time, there is an excess of power in the 
negatively polarized field at small k, indicating the presence 
of a large-scale field due to the mean-field dynamo, as also 
found in SB14. 

The blue curves at t = 500, in the top panel, show 
that the negatively polarized (large-scale) field (given by 
the solid line) near small k is starting to rise above the rest 
of the spectrum. Correspondingly we see from the bottom 
panel that the growth rate of the negatively polarized large- 
scale field component (solid curve) peaks at small k, whereas 
that of the positively polarized field (the dotted line) peaks 
at around k = 4. Furthermore, there is a decrease in A for 
the rest of the spectrum. This can also be noted from Fig¬ 
ure [S] where the A (k) at t = 500 for large k had decreased 
significantly as compared to when t = 200 or t = 300. Note 
that the peak in A (k) at k = 2 for the negatively polarized 
component, is consistent with the expectation that a fully 
helical a 2 mean-field dynamo grows fields at k = ki/2. 

Therefore, at these later times, there is no longer a grow¬ 
ing eigenfunction. Instead, the field is beginning to order 
itself both on large scales (due to the helical large-scale dy¬ 
namo) and also now on the forcing scale (due to the small- 
scale dynamo modified by the Lorentz force). Finally, at 
t = 900, shown as the red curves, the energy in the large- 
scale fields (E^ [ (k, t) given by the solid line) shows a peak at 
fci as seen in the top panel. In the bottom panel, we see that 
the growth rate for both polarizations are close to 0, but the 
negatively polarized large-scale field (solid curve) still shows 
a positive A at small k. Thus, at this time the small-scale field 
has saturated, but the large-scale field continues to grow in 
the presence of a saturated small-scale dynamo (albeit due 
to the loss of helicity through resistivity; see below). 

3.3 Evolution of the small-scale Held coherence 

As our work is focused on the evolution of the large-scale 
field in the presence of growing small-scale fields, it is also 



Figure 6 . Evolution of the integral scale L7) for Runs A, B, 
and C. In the inset, their respective normalized saturated final 
magnetic spectra are shown. 



Figure 7. Evolution of the integral scale /.) . for Runs D, E, 
and F. In the inset, their respective normalized saturated final 
magnetic spectra are shown. 

of interest to examine the evolution of coherence properties 
of the small-scale field. It is well known that the non-helically 
driven fluctuation dynamo generates fields whose power is 
concentrated on resistive scales in the kinematic stage. We 
have seen from our DNS (and from the work of SB14) that 
this continues to hold even when the turbulence is helical. 
In fact, in the kinematic stage, the spectrum is dominated 
by power which is concentrated at resistive scales due to the 
fluctuation dynamo. For a high-ReM system, if such a fea¬ 
ture persisted in saturation, the prevalence of a mean field 
would be questionable. Thus, it is important to investigate 
whether for a high-ReM system, there is a shift of magnetic 
energy from resistive scales to larger scales closer to the stir¬ 
ring scale on saturation. To address this question, we show 
in Figure [SJ the time evolution of the integral scale of the 
positively polarized spectrum (which is predominantly the 
small-scale field), for Run A. Similar results were obtained if 
we define the energy spectrum of the small-scale field to be 
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Figure 8 . Evolution of the integral scale, L. ., of the negatively 
polarized E^(k) characterizing the large-scale dynamo, is shown 
for Runs A, B, and C, corresponding to ReM = 330, 160, and 65, 
respectively. 



Figure 9. The ratio B/B rms is evolving with time, for Runs A, B, 
and C. In the inset panel, for Run A, the three curves of B rms (t), 
B (£) and y/kM^ for k = 4 are shown separately, where the latter 
two curves are scaled by a constant to make all curves overlap in 
the kinematic stage. 


with k > k{. This integral scale, defined here sep¬ 
arately for positively and negatively polarized components, 
is given by 


lL(!) 


f(2n/k)E±jk,t)dk 
f E+(k,t)dk 


(6) 


In the following, we are particularly interested in L[^ t (t), 
which characterizes the small-scale part of the field. In Fig- 
ure[Gj the integral scale for the fiducial high resolution Run A 
is shown as a solid black line. We also show the results from 
two lower resolution runs, Runs B and C, with PrM =0.1. 
In Figure0 we show the results from Runs D, E, and F, all 
of which have PrM = 10 . 

For the fiducial run, we see that the integral scale is 
roughly constant at L+ lt ~ 0.17 during the kinematic stage, 
reflecting the fact that the positively polarized field grows 
as an eigenfunction during this stage, with a small coher¬ 
ence scale. However, as the Lorentz force becomes impor¬ 
tant, begins to increase rapidly. This process begins 

at t ~ 400, which is also the time when the large k modes 
(k > 50) stop growing (see Figures Q] and 0. This rapid in¬ 
crease stops at t ~ 600 and L+ t ~ 0 . 8 , when the small-scale 
field modes with fc = fcf have largely saturated. There is 
then a subsequent slower rise of to ~ 0.9, as the large- 
scale field starts to dominate. Thus, there is considerable in¬ 
crease (by more than a factor 5) in the integral scale of the 
small-scale field from the kinematic to the saturated state. 
This factor is higher than that of ~ 3 seen for fluctuation 
dynamo in the purel y nonh elical case, albeit for PrM = 1 
(IBhat fc Subramanianl 120131 1. In the case of PrM = 10, as 
for Run D, ZA lt increases from a value of ~ 0.075 in the 
kinematic stage to ~ 0.6 at saturation, which is an increase 
by a factor of ~ 8 . 

Similar evolutions are also seen at lower ReM- However 
here L^ t is larger even in the kinematic state (reflecting the 
smaller resistive wavenumber for a lower Rcm). And Lf nt also 
saturates at a larger value for a lower ReM, reflecting the fact 
that the spectrum cuts off at smaller k for lower Rcm- This 
latter feature can be seen from the normalized spectra shown 


in the inset of Figure [ 6 ] Also note that, even though the 
saturated is slightly different for Runs A, B, and C, in 
all three cases, the peak power for the saturated small-scale 
dynamo is always at the forcing scale, fcf = 4. In fact, one 
can also compare the L^ nt obtained for our fiducial run with 
that expected for a small-scale field coherent on the scale of 
forcing. Suppose we modelled the spectrum of this field as 
E+ = M 0 {k/k{) 2 for k < k { and E+ = M 0 (fc/fc f )“ 5/3 for k > 
fcf, then one gets L[^ t = 0.6(27r/fcf) = 0.94, which compares 
reasonably well with that obtained in our highest resolution 
run. Thus, it appears that, on saturation, the small-scale 
held at the given forcing scale has become almost as coherent 
as possible. Note that, if the power in the small-scale fields 
were still at resistive scales, any peak in Em(k,t) at fc = 1 
would have made negligible contribution to B lms . Therefore, 
the above result provides some assurance that, even in high 
ReM systems, the large-scale held can indeed be signihcant 
and reveal itself on saturation. 

In the Figure [5] we also show the evolution of L ; “ t (for 
the large scale held) for the runs with PrM = 0.1. In the kine¬ 
matic regime, these curves are only slightly higher in ampli¬ 
tude compared to the corresponding lA it curves in Figure [fJ] 
thus indicating the slight excess in energy at large scales. 
Around the same time as Lf nt , the L;“ t curves start increas¬ 
ing to higher values. For Run A, Lr increases from 0.17 to 
a value of 2.5 which is about factor of 3 larger than the fi¬ 
nal Lf nt . The difference in the hnal values between the three 
curves for different ReM is because the large-scale helds are 
still growing due to resistive dissipation of small-scale heli- 
city. 


3.4 Growth of the large-scale field 

We define the strength of the large-scale held Bit), by inte¬ 
grating the energy spectrum Eyi{k, t ) between fc = 1 -2 and 

_2 

equating this to B /2. The ratio of strength of this large- 
scale held relative to the rms held B(t)/B Tms is shown in 
Figure [5] evolving from the kinematic stage to nonlinear sat- 
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Figure 10. The ratio B/B rms is evolving with time, for Runs D, 
E, and F. 


uration, for Runs A, B, and C. For the fiducial Run A, shown 
as a solid black line, B/ B rms ~ 0.04-0.05 in the kinematic 
stage. This is a factor 2 larger than B/B Ilns determined by 
SB14 (Fig. 6 ) on the basis of the mean power in horizontally 
averaged fields. However we find that in the kinematic stage, 
the ratio B/ B rma in fact decreases with ReM approximately 
as Re^ 3//4 , which agrees with the scaling found by SB 14 for 
the horizontally averaged mean field. 

The ratio B / B lma begins to increase rapidly once the 
small-scale field starts saturating, at t ~ 400 for Run A. 
This reflects the fact that, while the large k modes saturate, 
the k — 1, 2 modes which determine B, continue to grow 
due to large-scale dynamo action. This can also be seen ex¬ 
plicitly in the inset of Figure 0 where £> rms in solid black 
turns to saturate at around t = 400 while B in dashed blue 
continues to grow, again reflecting the continued efficiency 
of the mean (large-scale) field dynamo even when the fluc¬ 
tuation dynamo begins to saturate. Finally, after t = 600, 
this ratio enters a phase of slower growth; at this stage basi¬ 
cally the large-scale field becomes more and more dominant 
as the helicity of the small-scale field is lost due to resistive 
dissipation. Note that, although larger k modes begin to sat¬ 
urate by t = 400 as was seen earlier also in Figure [21 modes 
closer to the forcing scale (for example the mode in M 4 , also 
shown again in dashed red in the bottom panel of Figure [9]) 
are growing at the same rate as the modes in k = 1,2 be¬ 
tween t = 400 to 500 and they turn to saturate at nearly the 
same time, thus reinforcing the picture of a unified dynamo. 
But from t = 600 onwards, while the modes equal to k = 4 
and larger saturate completely, the k = 1,2 modes continue 
to grow on resistive time scales. By the end of our run we 
have B/B rma ~ 0.4, for the fiducial run. The lower ReM 
runs develop an even higher value of B/B rms ~ 0.5 and 0.6, 
for Rom = 160 and 65, respectively. This larger B/B lma for 
lower ReM arises because the resistive dissipation of helicity 
is more important in these cases. 

We show the evolution of the ratio B/B vma also in the 
case of PrM = 10 in Figure ITOl for Runs D, E, and F. Here 
again the ratio B/B Ilna increases dramatically from the kine¬ 
matic to the saturation regime. We can assess how the scal¬ 
ing of the ratio B/B Tlna vs. Re^ changes from the kinematic 


to the saturation regime. In the former case of PrM = 0.1, £ 
changes from a value of ~ —0.75 in the kinematic phase to 
~ —0.2 during saturation. I 11 the case of PrM = 10, the kine¬ 
matic phase is rather short in the current runs, but one can 
see that £ again tends to a small value of about — 0.2 upon 
saturation. The small residual ReM dependence of B/B Tma 
is also expected to disappear in the final saturated state of 
the a 2 dynamo, where one expects this ratio to tend towards 
(fcf/fci) 1 ^ 2 , provided the sm all scale helicity can be lost effec¬ 
tively jBrandenburdlioOlI l . In the current simulations as the 
helicity is only lost resistively, this has not yet occurred. All 
in all, we see that a significant large scale field can be gener¬ 
ated even in the presence of an active fluctuation dynamo. 
By the time they saturate the curves come closer together 
indicating that the mean field is dominant now. 


4 DISCUSSION 

Using high resolution DNS, we have shown that, for heli¬ 
cal turbulence at large ReM, the mean-field dynamo works 
efficiently to generate significant large-scale fields - even in 
the presence of a strong fluctuation dynamo. It appears that 
there is only one unified large/small-scale dynamo in such 
helical turbulence where initially fields on all scales grow 
together, and when the Lorentz force becomes important, 
successively larger scales saturate. 


4.1 Shape-invariant growth of the spectrum 

As in SB14, we find that in the kinematic stage the spec¬ 
trum grows as a shape-invariant eigenfunction of the helical 
dynamo, peaked at small scales (or large k ). There is clear 
evidence for a large-scale field even at this stage, seen as ex¬ 
cess power at small k in the negatively polarized component 
of the energy spectrum. However, the ratio of the strength of 
this large-scale field to the rms field decreases with increas¬ 
ing ReM- This is due to efficient fluctuation dynamo action, 
which amplifies power at small (nearly resistive) scales. The 
question then arises whether, in the presence of a fluctuation 
dynamo, the large-scale field can grow to a significant frac¬ 
tion of the rms field at least when the dynamo saturates. 


4.2 Scale-dependent saturation of the unified 
dynamo 

We show from the evolution of the spectra (Figures [T] and [2j 
that, as the field grows, small scales (large k modes, k > 4) 
saturate first, but the large-scale field (with k = 1 ) contin¬ 
ues to grow at about the same rate as the k = 4 mode, 
even when this happens. This can be seen by examining the 
wavelength-dependent growth rate of both the differently 
polarized components (Figure [4j and the total field (Fig¬ 
ure 0 . These A (k) start out as being independent of k in 
the kinematic stage, but progressively decrease to zero, first 
at large k and then at smaller and smaller k. This satu¬ 
ration behavior where small scales saturate first and then 
larger and larger scales saturate, remains qualitatively un¬ 
changed for PrM = 10 , even though the small scale dynamo 
is more efficient, as can be seen in Figure 0 
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4.3 Increase of small-scale field coherence 

At the end of our simulation, the spectra displays two peaks, 
one at the forcing wavenumber fcf, and the other at k = 1 . 
Therefore the back reaction due to the Lorentz force has 
enabled the small-scale field coherence to increase from small 
scales to the forcing scale, and at the same time allowed the 
large-scale field to develop. 

The first feature can also be seen from Figures [ 6 ] and [3 
where we show the evolution of the integral scale lA it of the 
positively polarized component (identified with the small- 
scale field). For our fiducial Run A (PrM = 0.1), we show 
that evolves from a value of ~ 0.17 in the kinematic 
stage to L+ t ~ 0.9 upon saturation, a significant fraction of 
the forcing scale ( 2n/k {). Also in the case of PrM = 10, as 
in Run D, increases by a factor of ~ 8 . In fact, through 
nonlinear saturation, the small-scale field has become as co¬ 
herent as possible for the given forcing scale. 


4.4 Significant growth of large scale field upon 
saturation 

The growth of the large-scale field to significant levels, even 
in the presence of the fluctuation dynamo, was also shown 
by considering the time evolution of B/B vms (see Figures [3 
and EB. This ratio, in the case of PrM = 0.1, starts from a 
small value of ~ 0.04 during the kinematic stage, but at the 
end of our run, we obtain a significant large-scale field with 
B /B lms ~ 0.4. A large increase in this ratio in seen for also 
Runs D, E, and F with PrM = 10. The growth of the ratio 
in Run A occurs in two stages: First between t ~ 400-600 
there is a rapid growth of B/B rms as the fluctuation dynamo 
saturates. It appears that the nonlinear ordering effects of 
the Lorentz force that saturate the small-scale fields, still 
allow growth of progressively larger scale fields, including 
scales larger than the forcing scale at k < 4. For t > 600, 
there is a slower growth of B/B Ilns presumably driven by the 
resistive dissipation of the small-scale helicity. It would be 
interesting to ask if, in very large ReM astrophysical systems, 
the effect of this resistive dissipation of small-scale helicity 
can al so be achieved by having instead magnetic helicity 
fluxes dBIackman fe Fieldll 2000 | : TKleeorin et al.ll 2000 f) . 


4.5 Quantum mechanical analogy 


It may be instructive to think in t erms of t he Kazant¬ 
sev model incorpo r ating helicity { Subramanfan| 1999|; 
iBoIdvrev et al.l 120051 : ISubramanian fe Brandenburel 12014) , 
where the dynamo problem is mapped to a quantum me¬ 
chanical potential problem, with growing small-scale dy¬ 
namo modes mapped to bound states in the potential and 
with helicity allowing for tunnelling to have enslaved large- 
scale field correlations with the same growth rate. The effect 
of the Lorentz force could be to make the potential well at 
the small scale k^ 1 to become shallower, allowing for only 
the marginally bound state to exist, while still having suf¬ 
ficient depth at the large scale fcf 1 , to allow the tunnelling 
free-particle states to grow. 

Such behavior is indeed obtained in a related real-space 
double-well potential problem, arising in non-axisymmetric 
galacti c dynamos, where the d ynamo is enhanced along a 
spiral (IChamandv et al.ll2013allbh . There the potential wells 


are near the galactic center and the corotation radius of the 
spiral. The fastest growing kinematic eigenfunction is largest 
in the central regions. But its tail is enhanced along the mag¬ 
netic spiral, near corotation radius. However, saturation of 
the field near the galactic center, still allows for the field to 
grow around corotation and become significant. From our 
work here, it appears that such a situation can also be ob¬ 
tained for a double-well potential in ‘scale’ or wavenumber 
space, when one incorporates nonlinear saturation effects. It 
would be of interest to demonstrate this also in a nonlinear 
version of the Kazan t sev m odel, perhaps generalizing the 
work of lSubramanianl (1 1994 : ferandenburg fe Subramanianl 
d 200 dl to include helicity loss. 


4.6 The very limited role of a effect growth rates 


In the early days of mean-field dynamo theory, computing 
linear growth rates was about the only way different dynamo 
modes could be characterized. In the late 1980s, their lim¬ 
ited usefulness became clear. Only the marginally excited 
case of zero growth rates remained truly useful. In partic¬ 
ular, the dominance of one mode over the other is entirely 
determined by nonlinearity, a nd not at all by differenc es in 
their kinematic growth rates IIBrandenburg et al.ll 19891 ). 

Linear growth rates have traditionally also been used 
to estimate the time it takes for the large-scale dynamo to 
reach saturation. The linear growth rate of the a 2 large-scale 
dynamo is expected to be much smaller than that of the 
small-scale dynamo. We have seen however that all modes 
grow together in the kinematic stage, and large-scale modes 
at k = 1-2 continue to grow at the same rapid rate even 
when the small-scale modes (k > 4) saturate. Note that it is 
the linear growth rate which has been important in discus- 
sio ns of the strength of the large-scale field in young galax¬ 


ies dKronberg et al.l 1992|; iBernet et al.ll2008l : Ijoshi fe Chandl 


l20ia T iFarnes et al.ll20l4 . Our present work does not really 

apply to galaxies, where shear is also important and the 
magnetic Prandtl number is large, but it highlights quite 
clearly that any estimate based on the value of a, or the 
value of |aV fl | 1;/2 in models with differential rotation Q. 
l|Beck et al .11 1 99)1 ) must b e irrelevant . This was in princi¬ 
ple already recognized by iBeck et alJ (1 1994h , who invoked 
a small-scale dynamo at early times to kick-start the large- 
scale dynamo at later times. 

In galaxies, the turnover time on the integral scale can 
be as short as 10 s or 10' years, but with fluid Reynolds num¬ 
bers well above 10 7 , the relevant e-folding time of dynamo 
growth will be shortened by a factor of Re 1 / 2 ~ 3 x 10 3 or 
more. Based on this argument, galactic dynamos may reach 
saturation first at the smallest eddy scales, on a time scale 
as short as several hundred years. A relevant limitation of 
reaching coherent large-scale fields comes only from the late 
saturation phase when magnetic helicity fluxes are expected 
to play an important role. 


5 CONCLUSIONS 

The results presented here support the idea that large-scale 
fields can be efficiently generated even in the presence of 
strongly growing fluctuations driven by the fluctuation dy¬ 
namo. Clearly, the growth of the larger scale field is aided 
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by the presence of helicity in the turbulence. But it is not 
as if there is an a 2 large-scale dynamo independent of the 
small-scale dynamo; as the growth rate of the k = 1 mode 
does not seem to change significantly right from the kine¬ 
matic to the nonlinear stage. Rather it appears that there 
is one unified dynamo, with all scales initially growing to¬ 
gether at one rate, and then the largest scales continuing 
to grow (aided by small-scale magnetic helicity loss) as the 
small-scale fields saturate. 

Several extensions of our model can be envisaged. Our 
dynamo is a homogeneous one, making catastrophic (resis¬ 
tive) quenching effects more pronounced. It would therefore 
be useful to extend our studies to inhomogeneous systems, 
for example when there is shear. In that case, the magnetic 
energy density of the mean magnetic field, B z , is no longer 
constant in space, which leads to a nonuniformity of the 
magnetic helicity flux divergence and can thereby allevi- 
ate catastrop h ic (premature) q uenching, as was shown by 
IHubbard fe Brandenburg (120121 1. Such models should there¬ 
fore be studied more thoroughly. However, as astrophys- 
ical systems are all confined in space with a corona and 
low-density material outside, it would be useful to address 
such systems in some fashion. It would be interesting to 
see whether this could remove the slow-down of the growth 
caused by total magnetic helicity conservation during the 
saturation phase. This would be particularly important in 
view of understanding the observed levels of coherent mag¬ 
netic fields in young galaxies. 
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